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Abstract 

The claim put forward in a recent paper by B. Jurado, K.-H. Schmidt and 
^ ' J. Benhiure that the transient effect of nuclear fission may be described simply 

(N , as a relaxation process in the upright oscillator around the potential minimum 

is refuted. Some critical remarks on the relevance of this effect in general are 



Tj- ' added. 

o ■ 



In the paper P it has been claimed that "a new, highly realistic analytical approxi- 
mation to the exact solution of the Fokker-Planck equation" has been presented. In 
this note we should like to raise some questions about the justifications of the ap- 

■ proximations used there, in particular with respect to its application to the decay of 
^ . a meta-stable state. Before we address details of the approach a few remarks of more 
^ ' general nature are in order. 

It has become customary to look at nuclear fission as a time dependent process. 
^ ■ As the current jh across the barrier shows a "transient behavior", simply because it 

■ takes some finite time before ji, reaches a quasi stationary value, one likes to interpret 
the result in terms of a time dependent decay width Tf{t), modifying in this way the 
one originally deduced by Kramers using essentially the same picture. This transient 
time seemingly implies a delay of fission during which light particles might be emitted 
in addition to those given by the conventional ratio r„/rf of the partial widths T^ 
for neutron evaporation to the stationary value Tf for fission, even if for the latter 
the Bohr- Wheeler expression Tbw is replaced by the smaller value Tk of Kramers' 
rate formula. In arguments in favor of such a procedure it is often claimed that 
the time dependence comes in only if fission is considered as a transport process 
underlying dissipative forces. In this one forgets that the transition state method 
is also based on "collective motion" which, in principle, like particle emission, is a 
time dependent event. It is only that for these processes one has become accustomed 
to apply widths calculated in a time m-dependent picture, in which, in addition, 
inherent averaging procedures are applied. Truth is that also Kramers' rate formula 
does not represent anything other than an inverse average decay time. This can 
directly be seen by exploiting the concept of the "mean first passage time" (MFPT). 
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For over-damped motion an analytic formula for the r^fpt can be derived from which, 
under the usual conditions for Kramers' rate formula, it follows that Fk = V'^mfpt, 
see e.g.pj. Interestingly enough, the value of the Tmfpt does not depend much on the 
initial position of the system, in clear distinction to the transient effect [B] • In |S] and 
|1] the concept of the MFPT has been applied to nuclear fission to examine, for the 
limit of over-damped motion, if more light particles may be emitted than given by 
the ratio F„/Fk. In [4] it has been demonstrated that Kramers' rate formula is valid 
only for simple potentials and under favorable conditions for the temperature. For 
potentials having some structure in addition to just one pronounced minimum and 
one barrier the fission lifetime was seen to be considerably longer than the tk = ^/Fk 
associated to Kramers' rate. This feature may already by inferred from the form 



the Tk takes on for over-damped motion. Any uncertainty in the product of the 
two stiffnesses at the minimum and the barrier, Ca and Cb, respectively, reflects 
itself in a corresponding error of their geometric mean, and hence in tk. Indeed, 
these stiffnesses are known at best in the immediate neighborhood of the extrema. 
Realistically, however, the potentials are hardly symmetric about these extrema. In 
cases that beyond the top of the barrier, for instance, the potential becomes wider 
this property may effectively imply a smaller |Cb| and, hence, a larger value of tk. 

We agree with the authors of [I] that the understanding of nuclear dissipation is 
of great importance, in particular its variation with shape and temperature. After 
all, these are perhaps the decisive features through which different models or theo- 
ries of nuclear transport can be distinguished It may perhaps be of interest to 
mention that, in addition to the papers cited in P, quite some work has been done, 
both experimentally jHl 13 as well as theoretically jH] , in which such questions have 
specifically been addressed. 

The main concern of ^ is that in previous analyses of experimental results un- 
certainties of a factor of two showed up in the so called reduced friction coefficient (3. 
There can be no question that the ultimate goal must be to improve our understand- 
ing about this problem, but it is questionable that the uncertainties and ambiguities 
of the method used in ^ imply progress. To begin with, one should not trace all 
problems back to just the one constant [3. Even discarding possible inaccuracies in 
the height of the barrier, which enters the decay rate in exponential fashion, there 
are crucial problems with the transport coefficients themselves. The /?, for instance, 
only stands for the ratio of friction 7 to inertia M. For obvious physical reasons, 
these two quantities must be expected to exhibit a totally different variation with 
temperature. Moreover, for a coordinate dependent inertia, Kramers' formula has to 
get an additional factor involving the square root of the ratio of the inertias at barrier 
and minimum [Hlini- For truly over-damped motion, on the other hand, any quantity 
which involves the inertia looses any meaning. Indeed, the latter does not appear in 
formula (P). 

Let us turn now to the more formal problems of |I] . The authors aim at delivering 
a simple way of calculating the time dependent pre-factor which supposedly relates 
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T{(t) to Kramers' stationary value Fk- The essential approximation is to calculate 
this pre-factor not from a global solution of Kramers' equation, which would properly 
account for the motion across the barrier, but from a solution of the same transport 
equation restricted to the upright oscillator by which the fission potential may be 
approximated in the neighborhood of the minimum. A moments reflection tells one 
that at the barrier, where the current jb is to be calculated, the height of the artificial 
potential in this region may easily exceed the barrier height several times. The very 
fact that in this region the stiffness of this auxiliary potential has the wrong sign is 
most crucial for the current, in particular at large times. Whereas for the inverted 
oscillator the jh{t) eventually turns into the stationary one already found by Kramers, 
the current for an upright oscillator tends to zero exponentially. In other words, 
replacing already in (pp|-6) (which is to say in eq.(6) of ref.p) the correct distribution 
Wn by the one for an upright oscillator, called W^^'^ in eq.(pP-9), leads to a vanishing 
denominator in this basic formula. 

To circumvent this problem some intermediate steps are performed to finally end 
up with formula ( lj-8) for which the W^'^^ of (pP-9) is to be inserted. One basic 
assumption for this is specified in eq.([T]-7). It implies that, for any time t and 
at the barrier top, the dependence of the distribution on coordinate and velocity is 
identical to the one at infinite time. For an oscillator this statement is easily seen to 
be incorrect, both for under-damped as well as for over-damped motion (for which it 
is claimed to be exact). Take the distribution given in (lT-9), namely^ 



which for an oscillator delivers the correct form for the density in coordinate space. 
For over-damped motion this statement is evident, for under-damped motion one first 
needs to integrate over velocity. Putting the form ^ into (P-7) the C{t) of (P-7) 
turns out to be 



which evidently is not only a function of time but varies with Xf,. As the upright and 
inverted oscillator turn into each other by analytic continuation (changing only the 
sign of the stiffness) |3] the proof just given also applies to the motion of a Gaussian 
across a barrier, if simulated by a parabola. 

For the cr^(t) needed for the W^^^{x = Xb, t) of eq.© a form is given in eq.(P-lO) 
which corresponds to zero initial width. On the other hand, the authors claim it to 
be more suitable to start from ground state fluctuations and they try to simulate 
this feature by introducing^ a time shift to- It is meant to represent the "time shift 

^It is properly normalized, also in the sense of eq.(l -5) if one only makes the common assump- 
tions that Xb is sufficiently far away from the minimum such that the tiny tail beyond Xb does not 
influence the normalization integral. 

^It remains unclear why the authors did not simple generalize the analytic form (pP-10) to one 
valid for any initial condition, such as ITJ, shown below for over-damped motion. 
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needed for the probability distribution to reach the width of the zero-point motion 
in deformation space", which is supposed to be "equal to the time that the average 
energy of the collective degree of freedom needs to reach the value (l/2)/ia;i associated 
to the zero-point motion". Obviously, the authors seem to understand to as a kind 
of relaxation time to the equilibrium of the oscillator, as represented by the ground 
state. It may be noted in passing that for a genuine quantum system any application 
of (PJ-10) is prohibited anyway as the distribution can never have zero width. In any 
case, it remains unclear why it should be the ground state and, hence, why there is 
no influence of the large intrinsic excitations which are produced in the first stage of 
the reaction; after all the authors work with a finite temperature. 

At this stage it may be worth while to remind the reader of some basic features of 
transport theory, which may help to clarify a few critical steps used in [T. To begin 
with, let us look how quantum features may be accounted for. As it stands, eq.(y[]- 
10) describes relaxation to the equilibrium specified by the equipartition theorem 
of classical mechanics, represented here by the pre-factor of the curly bracket. For a 
damped oscillator this may be generalized to represent quantum fluctuations correctly 
(see e.g. jS]). In equilibrium they are not given by those of the ground state h/2^uji 
used here, for instance in eq.([lj-12). In fact for over-damped motion just the opposite 
is true: there the correct quantum equilibrium is indeed given by the classical limit, 
see IS] for the oscillator and [TOj for the general case. 

Let us examine now the derivation of eq.(m-12), to = / {4:UiT) , meant to deter- 
mine the time lapse to- This equation is obtained by assuming a linear dependence 
between the cT^(t) and time t. This approximation is justified by arguing that the "in- 
fluence of the potential on the diffusion process" may be "neglected" as it "is anyhow 
small in the range of the zero-point motion". To see the catch in this argument let 
us simply write the correct equation for cr^(t), as it comes out of the Smoluchowski 
equation for the oscillator: 

^y{t) + 2-a\t) = 2D,^^. (4) 
at 7 

Here, C is the stiffness of the potential U{x), such that the latter may be written as 

U{x) = with C = fiul (5) 

and Dovd is the diffusion coefficient, which according to (jlj) is determined by the 
equilibrium fluctuation cr^ through 



I?ovd = -a2(t-oo)^-a,\ (6) 



Eq.(@]) implies that the "influence of the potential on the diffusion process" is given 
by the second term on the left which has the same size independent of the "range" of 
the coordinate. The solution of eq.(0]) for t > is given by 



0) -^eqjexp ( -- 



a\t) ={a\t = 0)- < ) exp ( t ) + <, (7) 
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showing that relaxation to the equihbrium value happens on the time scale Xovd = 
7/2C independent of the initial fluctuation. Of course, for t -C Tovd and zero initial 
fluctuations the cr^(t) becomes linear in t, but the reason why such a cr^(t) should be 
identifled as the ground state fluctuation of the undamped oscillator remains unclear. 
Put in the context mentioned above: it is unclear why such a value of cr^(t) = h/2fiuji 
should be relevant if reached by a process of strong damping. 

Next we turn to the value found for from eq.(rT]-12). It turns out so small that 
the introduction of this quantity and the associated fluctuation cannot explain why 
the r/(t) starts to become flnite only at about 0.7-10"^^ s. Indeed, for P = 2- lO^^s"^, 
hhJi = 1 MeV and T = 3 MeV one gets the very small number of to ~ 0.06 ■ 10~^^ 
s. The explanation given in [3] comes much closer: There, by simulating the whole 
fission process by a Langevin equation, it was demonstrated that such a shift is related 
to the relaxation of the initial distribution to the quasi-equilibrium in the minimum. 
For the numbers just used the relaxation time Tovd becomes Tovd — 0.36 ■ 10"^^ s — 
provided one makes use of the relation © . The value of Tovd becomes even very close 
to the ~ 0.7 ■ 10^^^ s at which in Fig.p-1 the T{{t) is seen to rise if the relation of 
frequency to stiffness is replaced by the m-correct one given in equation ([Tj-13) where 
the stiffness K is assumed to be only half the correct value given in (jSJ, in accord with 
the common definition used in text books not only on nuclear physics but on classical 
and quantum mechanics as well. It is true that, for the cases discussed in Fig. 
the motion is not really over-damped (for hui = 1), but the Tovd may nevertheless be 
taken fair estimate. 

As indicated before, it is left unclear why in the general case the system should 
start with a small fiuctuation. Indeed, formula (0-8) together with (p[]-9) implies any 
transient effect (of the type discussed here) to be absent if one chooses to start out of 
the quasi-equilibrium, which is to say for W^'^'^ ( t = 0) = W^°'^{x = Xb, t — > 00). 

For a bound system, like the upright oscillator, such an initial condition implies that 
the system stays in equilibrium for ever. For over-damped motion this may be seen 
from eq.(|7j) together with For a meta-stable situation like fission, on the other 
hand, the situation is different: then there will be a finite current outwards. Exactly 
this feature is not described correctly by formulas ([T]-8) and ([1-9). Please recall that 
in many cases an initial condition like that of the quasi-equilibrium specified before 
is not at all unrealistic. For sufficiently large fission barriers, as they are required for 
Kramers's rate formula anyhow, the system may well have enough time to reach such 
a stage around the first well before it decays by fission. 

Let us finally comment on the feature that for certain cases the present construc- 
tion seems to represent fairly well the numerically obtained global solutions of the 
transport equation for the full fission potential. In our opinion this feature should be 
considered accidental rather than supply a decent basis for trustworthy applications 
in future work of the approximations advertised in this paper. There are simply too 
many inconsistencies to warrant applicability to the general case. 



5 



References 

[1] B. Jurado, K.-H. Schmidt and J. Benlliure, Phys. Lett. B 553 (2003) 186 
[2] C.W. Gardiner, "Handbook of stochastic methods. Springer, 2002, Berhn 
[3] H. Hofmann and F.A. Ivanyuk, Phys. Rev. Lett. 90. 132701 
[4] H. Hofmann and A.G. Magner, to appear in PRC, see also nucl-th/0304022| 
[5] H. Hofmann, Phys. Rep. 284 (4&5) (1997) 137-380 

[6] D.J. Hofman, B.B. Back, L Dioszegi, CP. Montoya, S. Schadmand, R.Varma, 
and P. Paul, Phys. Rev.Let. 72, (1994) 470; see also: P. Paul and M. Thoen- 
nessen, Ann.Rev.Part.Nucl.Sci. 44, (1994) 65. 

[7] L Dioszegi, N.P. Shaw, I. Mazumdar, A. Hatzikoutelis and P. Paul, Phys. Rev. C 
61, (2000) 024613. 

[8] H. Hofmann, F.A. Ivanyuk, C. Rummel and S. Yamaji, Phys. Rev.C, 64 (2001) 
054316 

[9] V.M. Strutinsky, Phys. Lett. B 47, 121 (1973). 

[10] P. Petchukas, J. Ankerhold and H. Grabert, Annalen der Physik (Leipzig) (2000) 
1; J. Ankerhold, P. Petchukas and H. Grabert, Phys. Rev. Let. 87 (2001) 086802 



6 



